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Abstract 

We consider biological evolution as described within the Bak and Sneppen 
1993 model. We exhibit, at the self-organized critical state, a power-law sensi- 
tivity to the initial conditions, calculate the associated exponent, and relate it 
to the recently introduced nonextensive thermostatistics. The scenario which 
here emerges without tuning strongly reminds that of the tuned onset of chaos 
in say logistic-like onedimensional maps. We also calculate the dynamical 
exponent z. 
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There is nowadays a massive evidence of fractals and scale invariant phenomena in nature. 
They appear in an impressive variety of inanimate systems such as the geological (e.g., 
earthquakes) or climatic (e.g., atmospheric turbulence) ones, as well as of biological or 
living systems (e.g., biological evolution, cell growth, economic phenomena, among others). 
In most of the naturally occuring cases, no particular tuning is perceived. Per Bak and 
collaborators have advanced the hypothesis that, for many if not all the cases, this is so 
because the microscopic dynamics of the system makes it to spontaneously evolve towards 
a critical, scale-invariant, state. This is largely known today as self-organized criticality 
(SOC). Models have been formulated and experiments have been performed which profusely 
exhibit this interesting type of behaviour in sandpiles, ricepiles, earthquakes and others (see 
for instance and references therein). One such model is that introduced in 1993 by 
Bak and Sneppen P to paradigmatically describe biological species evolution. This is the 
model that we focuse herein. A variety of its properties are already known. Nevertheless, 
there is a crucial one that has never been addressed, namely the sensitivity to the initial 
conditions, which is known to be most relevant in nonlinear dynamical systems (quantities 
intensively studied such as Liapunov exponents, spread of damage, are in fact nothing but 
specific expressions of this concept). The study of this important property is the basic aim 
of the present work. 

Before describing our particular approach of this evolution model, let us introduce some 
preliminary notions by using, as a simple illustration, the following one-dimensional logistic- 
like map (see and references therein) 




This map recovers, for ^ = 2, the usual logistic map (in its centered representation). For 
fixed ( there is a critical value ac(C) such that , for a < adC), we observe a regular evolution 
(finite-cycle attractors), whereas, for a > ac{C), chaos becomes possible. Approaching ac{C) 
from below we can see the celebrated doubling-period road to chaos, with its successive 
bifurcations. Topological properties of the evolution (such as the nature of the successive 




1 - a\xt\\ (C > 1; < a < 2; t = 0, 1, 2, . . . 



(1) 
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attractors while varying a) do not depend on (, but metrical properties (such as the Feigen- 
baum's constants, characterizing geometrical rates of approach of the bifurcations leading 
to the chaotic critical point ac(C)) do depend on (Let us anticipate that the quantity 
we shall focuse herein, the entropic index q, belongs to this second class of properties). If 
we consider, at t = 0, two values of Xq which slightly differ by Ax(0) and follow their time 
evolution, we typically observe the following exponential behaviour for Ax{t) 

lim = exp[Ai t] (2) 

If Ai < (which is in fact the case for most values of a below adC)) we shall say that the 
system is strongly insensitive to the initial conditions. If Ai > (which is in fact the case 
for most values of a above ac(C)) we shall say that the system is strongly sensitive to the 
initial conditions. Finally, if Ai vanishes we shall speak of a marginal case. This is what 
happens, in particular, for the chaotic critical point. For this value of a, the sensitivity is not 
characterized by an exponential-law, but rather by a power-law. Eq. (2) can be generalized 
as follows @,|] 

We immediately see that the q = 1 case reproduces Eq. (2), whereas for g 7^ we have 

lim ~ [(1 - (t ^ 00) (4) 

Ax(o)^o Aa;(0) '^^ "J ^ ^ ^ ^ 

The g < 1 and the q > 1 cases respectively correspond to weakly sensitive and weakly 
insensitive to the initial conditions ||^. The coefficient Xg is a generalized Liapunov exponent, 
and satisfies Kg = Ag if Ag > and fCg = if Ag < (generalization, for arbitrary q, 
of the well known Pesin equality) , where the generalized Kolmogorov- Sinai entropy Kg is 
defined [|,|| analogously to the usual Kolmogorov-Sinai entropy {Ki herein). More precisely, 
in the same way Ki essentially is the increase per unit time of the Boltzmann-Gibbs-Shannon 
entropy Si = Pi^npi , Kg essentially is the increase per unit time of the generalized, 
nonextensive, entropic form [|l2 



S,^—^ (qen) (5) 

The nonextensivity of this form can be seen from the fact that if A and B are two independent 
systems (in the sense that the probabihties associated with A + B factorize into those of A 
and B), then 

S,iA + B)= S,iA) + Sg{B) + (1 - g) S^iA) S,iB) (6) 

We immediately verify that, since Sq is nonnegative, q = 1, q < 1 and g > 1 respectively cor- 
respond to the extensive, superextensive and subextensive cases. This generalized entropy 
(5) has generated a generalized thermostatistics [|l^,|ll|] (which recovers the usual, extensive, 
Boltzmann-Gibbs statistics as the q = 1 particular case), and has received applications in a 
variety of situations such as self-gravitating systems ||T^JT3|] , two-dimensional-like turbulence 
in pure-electron plasma ||T3| , [I^ , Levy-like [|l^] and correlated-like [1^ anomalous diffusions. 



solar neutrino problem JT^, peculiar velocity distribution of galaxy clusters [|T8l, cosmology 



||19|| , linear response theory ||2^, long-range fluid and magnetic systems [^, optimization 
techniques , among others (including of course the nonlinear maps (1) @J^). 



Now that we have introduced all the needed ingredients, we can close the illustration asso- 
ciated with Eq. (1) by mentionning that the (^-dependence of the entropic index q and of 
the fractal dimension df have been calculated |p at the chaotic critical point adO- It was 
exhibited that, while df varies from to 1 (or close to it), q increases from — oo to 1 (or 
close to it). Therefore, the Boltzmann-Gibbs limit g = 1 is attained when the attractor has 
an euclidean (nonfractal) dimension df = 1. But, when the system lives on a fractal {df < 1 
in the present case), nonextensive behaviour is revealed. 

Let us now return to the Bak and Sneppen evolution model. We shall exhibit that, at 
its self-organized critical state, weak sensitivity to the initial conditions (i.e., a power-law 
of the type indicated in Eq. (4)) occurs very similarly to the one just described for the 
chaotic critical point of the map (1). The model consists in a A^-site ring (linear chain 
with periodic boundary conditions); on each site {j = 1,2,. ..,N) we locate a real variable 



Bj {0 < Bj < 1, Vj) which corresponds to a "fitness barrier" separating two connected 
(first-neighboring) "species of hving organisms". We start by randomly and independently 
attributing the set of values {Bj}. At each successive elementary time step we identify the 
smallest Bj, and randomly change ("mutate") it as well as its two nearest neighbors. After 
some transient, a peculiar self-organized state emerges [0, rich in avalanches of all sizes and 
other scale invariant properties which makes the system to exhibit a variety of power laws. 
In the present work we have focused, as follows, the sensitivity to the initial conditions of 
the Bak and Sneppen model. Once SOC has been achieved, we consider that system as 
replica 1 {{Bj^^}) and create a replica 2 {{Bj^^}) by randomly choosing one of its sites, 
and exchanging the value associated with this site and that with the smallest barrier; we 
consider this moment as the collective time step t = 1 (we define a collective time step as A^ 
times the elementary time step, i.e., each site is going to be updated only once in average 
during a unit collective time step). From now on, we apply for both replicas the model 
rules of identifying the smallest barrier and updating that particular one and its two first- 
neighbors. We use (as in usual damage spreading techniques) the same random numbers for 
both replicas (hence, three different and independent random numbers are involved in the 
operation, since we change a particular barrier and its first-neighbors). 
We define now the Hamming distance between the two replicas as 

N 

N 



D{t)^W\Bf{t) - Bfit) (7) 



One such realization is shown in Fig. 1 for A^ = 1000. We then average A,, realizations (we 
have used typically N = 10^) and obtain {D){t): see Fig. 2, where a power-law is evident 
These results enable the determination of the slope 1/(1 — q) (see Eq. (4)), hence of q; we 
obtained q = —2.1 (to be compared, for instance, with the C = 2 logistic map value 0.24... 
IP). Finally, for fixed A^, if we denote by r the value of t at which the increasing regime 
crossovers onto the saturation regime (intersection, in Fig. 2, of two straight lines, namely 
those defined by the linearly increasing branch of the curve and the horizontal branch. The 
proportionality t{N) oc defines the dynamical exponent z ( and references therein). 



We obtained (see Fig.3) z = 1.56, to be compared, for instance, with 2.16 obtained for 
the square- lattice Ising ferromagnet. 

Since (g — 1) measures the degree of entropy nonextensivity (see Eq. (6)), it is an important 
index to be analyzed whenever discussing universality classes. Consistently, the determina- 
tion of q for other SOC models would be very welcome. Indeed, it will provide an insight 
on the fractal nature of the attractor towards which the system is spontaneously driven. 
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FIGURES 

FIG. 1. Time evolution of the damage associated with one realization of a typical system with 
N = 1000 . 

FIG. 2. Average of Nr realizations such as those of figure 1 for three different sizes: N = 1000 
(top), N = 500 (middle) and N = 250 (bottom). The slope 1/(1 - q) equals 0.32 ± 0.01. 



FIG. 3. Log-log plot of r versus N for the three curves of figure 2 . 
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